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_i ' Over the last years, Effective Field Theory (EFT) methods have become the standard tool to deal with strong 
interactions in the non-perturbative regime. In the form of Chiral Perturbation Theory (ChPT), they have been 
used with some success in both the mesonic and single nucleon sectors. The EFT of nuclear forces based on a chiral 
expansion was originally suggested by Weinberg pj, and since then a lot of work has been devoted to gain a better 
understanding of the two- nucleon interaction at low and intermediate energies While at the beginning, these 

\ studies did not aim at substituting the highly successful realistic potentials built from meson exchanges (Bonn-Jiilich, 
Nijmegen, Argonne, • • • potentials), the importance of uncovering such an EFT cannot be ignored, since this theory 
will allow for rigorous calculations of both elastic and inelastic processes in systems with two or more nucleons, in 
a framework consistent with the Standard Model of strong and electroweak interactions. Furthermore, the latter 
EFT works provide an accurate description of nucleon-nucleon (NN) phase-shifts for several partial waves and in a 
wide range of energies fl6| . Weinberg's original proposal was to determine the NN potentials using the organizational 
principles of ChPT and then to insert these potentials into the Lippman-Schwinger Equation (LSE) to solve for 
NN scattering amplitudes. Though such a scheme has proved to be successful 0], 0, 0, it suffers from formal 
inconsistencies, in particular, divergences that arise at a given order in the chiral expansion cannot be absorbed by 

i-^ | terms of the same order in the Lagrangian The trouble arises, since Weinberg's power counting assumes that 

the coefficients of the higher order operators were set by a typical Quantum Chromo Dynamics (QCD) scale, but in 
practice it turns out that there appear explicit M factors (being M the nucleon mass) in the needed counter-terms 
to renormalize the amplitudes, which formally decrease the degree of the operator in the Weinberg's counting. This 
was partially solved by Kaplan, Savage and Wise (KSW) who introduced a new scheme where pions are treated 
perturbatively Q- However, it turns out that the KSW expansion converges slowly in the 1 5o ( 2I+1 Lj) channel and 
does not converge at all in the 3 Si— 3 D± C0U Pl e( i channels |(J, |ld| . Very recently, a so called expansion about the 
chiral limit has been proposed in Ref. [l3| . which seems to be equivalent to KSW power counting in the 1 S'o channel 
^ . and to Weinberg power counting in the 3 Si — 3 D\ coupled channels. 

In this work we adopt the following point of view. For a scattering process involving external momenta < q, one 
should only consider a Lagrangian/potential which explicitly includes light degrees of freedom for which m < q. The 
effects of heavy virtual particles appear as an infinite number of non-renormalizable operators suppressed by some 
mass scale relevant to the degrees of freedom excluded from the theory. At low-intermediate momenta, below twice 
the pion mass, it should be enough to include explicitly one-pion-exchange, and simulate the rest of the physical 
contributions by a tower on non-renormalizable contact interactions, organized as a derivative expansion. Given a 
potential, which includes contact terms up to some order, the LSE performs an infinite sum of diagrams, which in 
general would require to be renormalized. After renormalization the coefficients of the contact terms would contain, 
in general, contributions from all orders in m T (pion mass). 

Most of the formal problems appearing within the Weinberg's scheme are linked to the use of a somehow restrictive 
scheme to renormalize the LSE Such a restricted conception of the renormalization of the LSE would also 

lead to unexpected consequences in other scenarios, as for example 7T7t scattering in the p— channel Among 
the abundant literature on the subject, the framework presented in Ref. |8| constitutes a first step towards the 
renormalization scheme proposed in this work. There, a subtracted LSE is derived, and the numerical values of the 
fitted Low Energies Constants (LEC's) are not translated to values for any Ultraviolet (UV) cutoff. 

The aim of this paper is to present a renormalization scheme of the LSE for NN scattering by a combination of known 
long-range and unknown short-range potentials. For simplicity, we restrict the discussion to the X S§ channel, where 



2 



we show that the proposed scheme describes phase-shifts up to Center of Mass (CM) nucleon momenta of the order of 
260 MeV. This renormalization scheme is based on previous findings on tttt 0] and meson-baryon scattering (l9j-|2l^. 
it leads to renormalized amplitudes which fulfill exact elastic unitarity, and it is easily generalized to study higher 
waves and/or number of derivatives in the contact part of the interaction. 



II. EFFECTIVE POTENTIALS AND THE LSE 



After projecting into the 1 S'o partial wave, the NN LSE for a CM nucleon kinetic energy E, reads: 

T(E;p,p) = V(p,p) + J Q dqq 2mE _ q2+ie 0) 

with to = 469.46 MeV the NN reduced mass, for which we take that of the neutron-proton system, p and p' are the 
initial and final relative momenta of the two nucleons, and V(p,p') the 1 Sq NN potential. The normalization of the 
scattering amplitude T(E;p,p') is such that on the mass shell, p = p 1 = \llmE = fc, it is related to the phase shifts, 
S, by: 

9 P 2i<5(fc) _ -, 

T(k) = — — — — - (2) 

7T 2lfc 

The potential consists of the one-pion-exchange potential, V n , plus contact terms, V s , with up to two derivatives Q: 

V(p,p') = V s (p,p') + V n (p,p') 
V s (p,p') = g a + gi(p 2 +p' 2 ) 

T/ f n 2ma^ f +1 dx ma^ 1 p 2 _ + to 2 g\m\ 

V*{p,p) = / 2 /2 — — ■ 2" = -log-g— — with a^ = -^— (3) 

7f 7-i p +p —2pp'x + m% tt pp p++m% ±vTrj£ 



with p± = p ± p', m K = 138 MeV and / w = 93 MeV the pion mass and weak decay constant and finally gA = 1-25 
the axial coupling constant. Note that the contact term includes the 5 3 (r ) contribution from one pion exchange, 
and it also does the leading and sub-leading contributions in the derivative expansion of all shorter distance effects, 
such as Two Pion Exchange (TPE) 1 , intermediate A's, u> exchange, etc • • •. Such a procedure suffers from some 
limitations, for instance, since the TPE potential has not been explicitly included one cannot relate NN scattering 
to other processes as pion-nucleon, pion-deuterium, etc . . .scattering. The LEC's go and g\ are not fixed by chiral 
symmetry and have to be determined by a fit to the phase shifts, as we will see below. Scattering by short-range 
interactions in the presence of a known long-range potential,^, can be treated by Distorted Wave Theory (DWT) 15]. 
We write the full scattering matrix as 2 

T = T 7r + (I + T n G )f s (I + G T 7r ) (4) 

with I the identity, Go the free nucleon Green function and T w the scattering matrix for V v alone. Besides, T s describes 
the scattering between distorted waves of . It satisfies the LSE 

t, = V s + VsG^Ts (5) 

where is the nucleon Green function in the presence of V n , i.e., G^ = Go + GqT^Gq. To solve the LSE of Eq. 
the full off-shell scattering matrix T^(E;p,p') is required, which is determined by solving the LSE of Eq. with 
the obvious substitution of V — ► V w . This equation is finite (it corresponds to the scattering by the usual Yukawa 
force studied in most books of Quantum Mechanics) and does not require to be renormalized, being then possible 



1 Note that, the full structure of the logarithmic terms, stemming from the pion loops implicit in the TPE contribution, cannot be entirely 
accounted for the contact terms. Indeed, for CM transferred momenta above 2m,, TPE contributions would have to be treated in the 
same footing as the one-pion exchange ones. 

2 Since the angular integrals have been already performed when projecting into the 1 So wave and taking into account that do not lead 
to UV divergences, in what follows we will assume that the measure in momentum space and the free nucleon propagator are given by 

So"" ^ G 0( E 'P'P') = S{P ^ P ] 2mS-p2+ie ' res P ectivel y- 



3 



a numerical evaluation. We have obtained T„ by discretizing the momentum space and using the inverse matrix 
algorithm. Some results are shown in the top panel of Fig. ^ To obtain T s (E;p,p'), we solve 

rf (J? f \ v( „, f + °° , 2 Vs(p,q)f s (E;q,p') , q 2 q' 2 'V s (p, q)T n (E; q, q')f s (E; q> ,p>) 
Ts{E;p,p) = V s {p,p)+ / dqq -= 5— h / dqdq -j- , ; . Wf2 ,, : . , (6) 



k 2 -q 2 + ie J Q (k 2 - q 2 + ie) (k 2 - q' 2 + ie) 

The above equation can be reduced to a linear algebraic system of equations, which solution is straightforward, 

f a (E,p,p') = a + 0(p 2 +p' 2 )+ 1P 2 p' 2 , (7) 

with the energy dependent functions 

go + g 2 K 4 3i - g\K 2 g 2 K 

a = t , 0= t , 7= — t — 

A A A 

A = I - g K Q - 2 9l K 2 + g 2 (K 2 - K K 4 ) . (8) 
The UV divergent integrals K n {E) = I n {E) + J n (E) are defined by: 

f°° a n+2 

W) = / a , . (9) 

J k z - q 2 + ie 

T (m f 00 dqdq'q a q ' b T«{E;q iq ') 

n[ ' ~ Jo (k 2 -q 2 + ie) (k 2 -q' 2 + ie) (iUj 

and the pair of integers (a, b) take the values (2, 2), (4, 2) and (4,4) for the n = 0, 2 and 4 cases respectively. 

We use the Dimensional Regularization (DR) scheme, since it preserves chiral and gauge symmetry and Galilean 
invariance, which makes the integrals relatively simple to evaluate. DR discards all power-law divergences of the type 
J dqq n , which makes finite all the I n integrals define above, 

7 B (£) = Fxi(— j (11) 

To simplify the J n integrals it is useful to realize that in DR the linearly UV divergent integral / + °° dqq 2 V R (q,p) is 
finite (= 2ma 7r m 7r ) and independent of p Making use of the LSE satisfied by to get expressions involving the 
above integral, the J 2 and J4 integrals can be related to the Jo one, 

= (E- m^)J n . 2 {E) + ^-m^k n - x (12) 

for n = 2,4. Besides, Jq(E) is logarithmically divergent and it only requires one subtraction, i.e., Jq(E) = Jo(E) — Jo(0) 
is finite and can be numerically evaluated. Plugging Eq. (J7J) into Eq. J3J one gets for the on-shell scattering amplitude 
T(E;k,k) = T{k) : 

T(k) = T 7T {k)+f s {k) + 2{{a + 0k 2 )L o + (l3 + 1 k 2 )L 2 }+aL 2 + 20L o L 2 + 1 L 2 (13) 

with T v (k) and T s (k) the long and short range on-shell matrices and 

[°° dqq n + 2 TAE;k,q) 
Ln{k) = Jo (k 2 -q 2 + ie) (W) 

The integral Lq is finite and the UV divergent integral L 2 in DR becomes finite, i.e., L 2 (k) = 2mx 
{(E — a 7r m^)Lo(fc) — a^m^. With all above results one gets 

T{k) = T^k) + (1 + L (k)) 2 f s {E- k, k) = T n {k) + ^^fr^ 77^ (15) 

V s (k, k) + mk/2 - J (E) - J o (0) 

with k = k 2 — 2ma 7r m 7r . Elastic unitarity is exactly fulfilled thanks to a Watson's type relation satisfied by Lo(fc) 3 , 

L (k) = (l (k)-mk/2)T 7r (k), l Q eH (16) 



3 This relation is easily obtained from the Optical Theorem satisfied by T n , i.e., 2ImT 7r (_E;pi,p2) = —T^kT lv {E;p\,k)T*{E;p2,k). 
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and that 

ImJo(fe) = - {(7rfc/2) 2 - Z (^o + 2RcT- 1 (fc))} lmT„(k). (17) 
Indeed, the on shell scattering matrix, T, can be re-written in a form where unitarity is manifest, 

V ' n W ReT^\k)+g{k) V ' 

ff(fe) - V-\k, k) - ReJ (fc) + l2 J 5°!(k) (19) 

Above threshold, the function g(k) is real, and <5„-(fc) are the phase shifts deduced from T„(k). The functions lo(k) 
and ReJo(fc) and the phase shifts are plotted in the top panel of Fig.Q] Note that in the complex E plane, g(k) 
does not have right-hand-cut, but it does have a left-hand-cut due to pion exchange. Before going ahead we should 
study how to get rid of the UV divergences. 



III. RENORMALIZATION SCHEME. 



The DR scheme has led to a drastic reduction of UV divergences and thus we are just left with the logarithmic 
divergent integral 4 : Jo (0), which should be absorbed into a redefinition of the coupling constants of the potential, 
provided one included in the potential all terms consistent with the symmetry principles. However, the amplitude 
cannot be made finite by simply redefining the couplings go and gi . This problem arises because we have not included 
operators with more than two derivatives in the potential, which are needed as counter-terms to render multiple 
insertions of the two derivative operator finite. Indeed, the divergent constant Jo(0) appears in the function g(k) 
defined above, and it seems natural to define a renormalized contact potential Vf 1 (k, k) — V ! T 1 (k, k) — Jq(0) which 
leads to an infinite series of powers of k 2 for V R . Since the LSE performs a non-perturbative resummation of diagrams, 
in principle one should also add a series of counter-terms to the bare amplitude such that the sum of both becomes 
finite. At each order in the k 2 (or derivative) expansion, the divergent part of the counter-term series is completely 
determined. However, the finite part remains arbitrary. It means that the coefficients of each of the terms in the k 2 
series, implicit in the definition of V S R , become undetermined and have to be fitted to data or, if possible, evaluated 
in QCD. Thus, the scenario might look like pessimistic, and because of the appearance of divergences, g(k) turns out 
to be a completely undetermined function. This situation has some resemblances to ChPT in the meson-meson and 
meson-baryon sectors. For simplicity, let us pay attention to elastic SU (2) irir scattering. It is well known |23| that 
the divergences arising at one loop, C(p 4 ), cannot be absorbed into a redefinition of the leading terms 5 , 0(p 2 ), of 
the Lagrangian. New counter-terms, U, higher in the chiral expansion are needed to absorb the divergences, which 
finite parts h, remains undetermined, cannot be fixed by chiral symmetry and have to be fitted to data. These 
LEC's are fundamental parameters of the EFT, which contain the contribution at low energies of higher degrees of 
freedom, which have been integrated out. The evaluation of the divergent parts of the amplitudes with an UV cut-off 
does not necessarily provide a reasonable estimate for them 6 . To be predictive, one can adopt here a renormalization 
scheme, similar to that used in Refs. |l8|-|2dj to renormalize the Bethe Salpeter equation for meson-meson and meson- 
baryon scattering, which reduces the enormous freedom discussed above. Since for small momenta, higher derivative 
operators should have a tiny influence, we choose to reduce all this proliferation of LEC's, by imposing relations, 
among all LEC's associated to counter-terms with a number of derivatives higher than the terms included in the 
potential, in such a way that the renormalized amplitude can be cast, again, as in Eq. I|18H and therefore it fulfills 
elastic unitarity. This amounts in practice, to interpret the previously divergent quantity Jq(0) as a renormalized 



4 



r. 



The divergent part of Jo(0) is given by the integral 



I — d / IT 7TL \ 

dqdpV-n(q,p) = —2ma- n I — ( arctan J . (20) 

o Jo 1 V2 1 1 

Power divergences can be absorbed into a redefinition of the leading terms, but logarithmic ones cannot, and require the inclusion of 
higher order structures. 

If an UV cut-off A is employed, the C(p 4 ) contributions in all isospin-angular momentum elastic tttt scattering channels will be determined 
just by one parameter, A, while at this order there are four independent parameters, ^i ,2,3,4 , which, for instance, incorporate t he e ffect 
of the p and other resonances in the amplitudes (see for instance discussion in Section 3 (pages 63—70) of second entry of Ref . Il8l ) . 
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undetermined parameter. After having renormalized, we add a superscript R to differentiate between the previously 
divergent, Jo(0), and now finite quantity 7 , J R (Q). This parameter and therefore the renormalized amplitude can be 
expressed in terms of physical (measurable) magnitudes. The estimate given by means of a reasonable UV cut-off for 
the numerical value of Jq (0) might not be good. For instance, Eqs. (A15) of second entry of Ref. [18( illustrate this 
point, if an UV cutoff is employed, the divergent integrals appearing there will be independent of the I J channel, and 
to get a reasonably value for unrealistic scales or cut-offs of the order of 300 GeV will be needed. This is due 

to the special nature of the p-resonance, and that is not the case, for instance, for the s— waves. 

Thus, we are proposing an expansion of the short range part of the potential in powers of fc 2 , i.e., delta function 
and its derivatives in coordinate space, and for a given potential, to compute T to all orders in the fc 2 expansion to 
restore exact unitarity. Thus, the contact terms effectively account for degrees of freedom higher than the pion and 
one should expect the scheme to work up to the first threshold, likely two pion production {k^AB ~ 2m w ), which 
is around the normal nuclear matter Fermi momentum. A word of caution must be said here; this expansion is not 
equivalent to a chiral expansion in m n , since the coefficients of the contact interaction would contain contributions 
from all orders in m T , as it was firstly pointed out in Ref. |^|. 



IV. RESULTS AND CONCLUDING REMARKS 

After the above discussion, it is clear we have three undetermined LEC's: go, g\ and Jq (0), which we obtain from 
a x 2 — fit to Nijmegen phase-shifts [24|. To perform the fits, we assume errors on the phase shifts, S (in degrees), 
given by 0.5 + abs(<5) x max (0.1, E^a^MeV]) /100, where is the kinetic energy in the laboratory frame. We fit 
data from threshold up to a certain value, fc max , of the CM nucleon momentum. We determine fc max , by studying the 
dependence of x 2 /dof on it, and fix it in the region of 260 MeV, since the inclusion of simply a new datum of higher 
energy would double the value of x 2 /dof . Thus, we have fitted 62 phase shift values for which E^ a \, has varied from 
0.01 to 142 MeV, and results are shown (solid line) in the bottom panel of Fig.Q] Best fit parameters are 

J R (Q) 

g a m = -0.276(6), gi m 3 = 0.347(14), -2-^2 = -3.21(8), (21) 

m 

with x 2 /dof = 0.016. In the above equation, statistical errors have been shown in brackets, affect to the last digit of the 



parameters and are given by the square root of the diagonal elements of the covariance matrix, v. 



i j 



\2db k db, J 



with bi any of the three fitted parameters. It is worth to mention that from the above matrix, one also learns that 
the pairs go,gi and go, Jrj^O) are highly anticorrelated, with correlation coefficients (r^ = Vij / \ZvuVjj) smaller than 
—0.99. This means that effectively there is just one independent parameter. We will come back to this point below. 

If the derivative terms of the contact interaction are set to zero, g\ = 0, the UV divergence can be absorbed in go 4] 
and thus one is just left with one parameter, which can be fitted to data. Yet, though g\ ^ 0, one can arbitrarily set 
the renormalized coefficients of the higher order terms to zero, i.e., take (0) = 0, and fit the non-zero LEC's, go 
and gx, to data. Both procedures are also shown in Fig.^ the first one leads to the scheme developed in Ref. @ and 
the second one was studied for the very first time in Ref. As seen in the figure, these two schemes do not work 
for momenta higher than 10 and 40 MeV respectively, while the solid line prov ides a good description of data up to 
momenta of the order of 260 MeV. Indeed, from Eq. I|21(l we have \/|.9o/2.9i| ~ 295 MeV, while such ratio takes a 
value around of 70 MeV when J(f (0) is set to zero. 

The scattering length, a, and the effective range, ro, are defined from the effective range expansion, 

T-\k) = ~ (-- + \r Q k 2 + ik) (22) 



2 V a 2 

Fitting the inverse of our amplitude, up to k = 48.5 MeV, to the above formula, we get a = —23.65 fm and ro = 2.63 
fm, in reasonable agreement with the Nijmegen results ,24] (—23.7 fm and 2.73 fm, respectively) also obtained in this 
way. 



7 Thus, in the calculation presented in this work, we relate the finite parts of the k 6 , k s , etc., contributions to that of the fc 4 one 
(note that, the fc 4 contribution shows up at the one-loop level induced by contact terms quadratic in momenta: p 2 and p' 2 ), which is 
determined by J^(0). The relation is such that elastic unitarity is restored, and thus this scheme differs from those where the higher 
order terms are set to zero. Indeed, setting to zero the higher order terms is as arbitrary as setting them to any other value. At 
next order in this derivative expansion, it is to say when terms of the type p 2 p' 2 , p 4 or p' 4 are included explicitly in the potential, 
the LSE would generate terms of order fc 6 and higher in the amplitude. Thus, the expansion based on this renormalization scheme is 
systematically improvable. 
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FIG. 1: Top panel: Several quantities extracted from T^: phase shifts and functions lo(k) and ReJo(fc) = ReJo(fc) — Jo(0) defined in 
Eqs. 1121 and 1161 . respectively. From the energy dependence of &^{k), we deduce the effective range parameters (see Eq. 1221 1 a-n- = —0.88 
fm and ro^ = 12.38 fm. Bottom panel: 1 So np phase shifts in degrees plotted versus CM momentum (log scale). Solid line stands for the 
results of this work (Eq. 1211 '). while other two approaches with J^(0) = and J^(0) = gi = respectively, are also shown. Data are 
taken from the phase-shift analysis of Ref. 12411 . and data for momenta above the vertical line (260 MeV) have not been included in the fit. 



To finish, we would like just to summarize the results obtained in this work. We advocate for an expansion of 
the short range part of the potential in powers of k 2 and computing T to all orders in the k 2 expansion to restore 
exact elastic unitarity. DWT techniques and DR have been used to solve the 1 S$ NN LSE with a potential, which 
consists of pion exchange and contact terms with up to two derivatives. The procedure of solving the LSE is quite 
simple and an explicit expression, where exact elastic unitarity is manifest, has been given as well (see Eq. O). 
A particular renormalization scheme has also been discussed. It is based on previous findings on meson-meson and 
meson-baryon systems and its main ingredient is to realize that an EFT is not a renormalizable field theory in the sense 
of QCD, i.e., with a finite number of counter-terms. Thus, to keep finite the amplitude, obtained after performing 
the non-perturbative resummation implicit in the LSE, would require, in general, the addition of an infinite set of 
counter-terms in the short distance potential (k 2 series). The finite parts of the coefficients of the series would remain 
undetermined (LEC's) and encode the contribution of higher degrees of freedom, not explicitly included. An UV 
cut-off can effectively account for this freedom only up to some momentum scale, which will depend on the terms 
explicitly included in the potential 8 . For a contact potential including up to k 2 terms, we find a good description of 



The EFT only becomes cut-off independent when all counter-terms compatible with the underlying symmetry are included in the 
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phase shifts up to CM NN momenta of the order of 260 MeV with the inclusion of just one additional parameter, 

J ° R{0) - 

In Ref. [7j, an UV cut-off is used to regularize the amplitudes, and with reasonable values of it, in the range 
0.6 — 1 GeV, a good description of data is found. Note, however, that this does not have always to be the case 
and depends, as we mentioned above, of the physical system and of the order of the expansion included in the 
potential. Thus, to describe the p— resonance in irir scattering, one is left to deal with UV cut-off's of the order of 300 
GeV ^3;E3- Unrealistic values of the UV cutoff are also needed to account for the iV(1650)— resonance in meson 
baryon scattering 0],[2(]], but this is not the case, for instance, for the a and /o(980) resonances in mr scattering [2l| 
or the A(1405) resonance in meson baryon scattering p2l |. 

Besides, we find extremely big anticorrelations between J(f(0) and the coefficients, go and <?i, of the iterated short 
distance potential and among the two later ones, as well. This indicates that the higher derivative operators in the 
EFT, generated in this scheme by the inclusion of J^(0), are actually highly correlated. This might support the idea 
that though, the higher derivative operators are controlled by a scale that diverges as \a\ — > oo, thanks to these high 
correlations, the effects that diverge with a cancel, as it was pointed out in Ref. Q. 

To include more derivatives is trivial and would result, besides the left hand cut due to pion exchange and accounted 
for by means of T n , in a higher order pade approximant for the function g(k) in Eq. (|18|l . The extension of the procedure 
to higher partial waves is also straightforward although cumbersome and will be presented elsewhere |25) . 
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